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Various aspects of the high-energy emission from relativistic jets associated with compact 
astrophysical systems are reviewed. The main leptonic and hadronic processes responsi- 
ble for the production of high-energy 7-rays, very-high energy neutrinos and ultra-high 
energy cosmic rays are discussed. Relations between the 77 pair production and photome- 
son production opacities are derived, and their consequences for the relative emission of 
7-rays and neutrinos are examined. The scaling of the size and location of the various 
emission zones and other quantities with black hole mass and dimensionless luminosity 
is elucidated. The results are applied to individual classes of objects, including blazars, 
microquasars and gamma-ray bursts. It is concluded that if baryons are present in the jet 
at sufficient quantities, then under optimal conditions most systems exhibiting relativis- 
tic jets may be detectable by upcoming neutrino telescopes. An exception is the class 
of TeV blazars, for which 7-ray observations imply neutrino yields well below detection 
limit. 
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1. Introduction 

It is widely believed that the high-energy emission observed in several classes of 
Galactic and extragalactic sources, e.g., blazars, microquasars, gamma-ray bursts 
(GRBs), is associated with coUimated, relativistic outflows. This view is strongly 
supported by some recent exciting discoveries, most notably (i) the discovery that 
many of the compact extragalactic radio sources are extremely luminous, hard 7-ray 
sources, (ii) detection of superluminal motions and very-high energy 7-ray emission 
in several Galactic X-ray transients, (iii) indirect evidence for jets in GRBs, and 
(iv) detection of a mildly relativistic outflow from a magnetar following a giant flare 
(GF). An overview of the observations is given in Sj2] below. 

A common view is that relativistic jets associated with compact astrophysical 
systems are powered by a magnetized accretion disk and a spinning black hole and 
coUimated by magnetic fields and/or the medium surrounding the jet. The outflows 
associated with highly super-Eddington sources are presumably driven by violent 
events, specifically collapse of a massive star or coalescence of compact objects in 
the case of GRBs, and magnetic field annihilation in the case of magnetars. Unfortu- 
nately, the mechanism responsible for the formation, acceleration and coUimation of 
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the jets is poorly understood at present. Even the composition of the jet is unknown 
in most sources. The high-energy emission exhibited by those systems is beheved 
to be produced close to the inner engine, in regions where dissipation of the bulk 
(kinetic plus magnetic) outflow energy occurs and, therefore, provides an important 
probe of the inner jet zone. As discussed below, some constraints on jet parameters 
can be imposed from the observations. Further progress in our understanding of 
the physics of relativistic jets and their engines is expected in the coming few years 
with the advance of observational techniques. A new generation of experiments just 
started operating or will become operative soon: Space-based (e.g. GLAST, INTE- 
GRAL, AGILE), Cerenkov (e.g. HESS, MAGIC, VERITAS), and air-shower (e.g. 
MILAGRO) 7-ray detectors will probe with high sensitivity the energy range of 
10 MeV to a few TeV (see, e.g., Ref. [T]for reviews). The operating BAIKAL and 
AMANDA neutrino telescopes, and the cubic-km scale telescopes under construc- 
tion (IceCube, ANTARES, NESTOR, NEMO; see, e.g. Ref.Elfor review), wiU open 
a new window onto the Universe. Besides providing an important probe of the in- 
nermost regions of compact astrophysical systems, these experiments can also be 
exploited to test new physics. Finally, new ultra-high-energy cosmic-ray detectors, 
e.g. the HiRes ^ and the hybrid Auger detectors will provide cosmic-ray data of 
unprecedented quality and quantity. 

1.1. The Basic Model 

The system under consideration is depicted in Fig. [T] The main ingredients relevant 
to the discussion that follows are indicated. In blazars and microquasars radiation 
emerging from the disk (represented by the dashed lines in Fig. [T]) may be an 
important source of seed photons in the jet. In GRBs the disk is highly opaque to 
electromagnetic radiation and cools through emission of MeV neutrinos. Moreover, 
the jet must first break out of the envelope of the progenitor star, driving a strong 
shock. The break-out shock may give rise to additional emission of 7-rays and 
neutrinos at early times, as will be described below. In microquasars a companion 
star may be present, as illustrated in Fig. [2] Its radiation and stellar wind can 
interact with the jet and may contribute additional opacities. 

2. Overview of Observations 
2.1. Blazars 

Many compact extragalactic radio sources, called blazars, have been found by the 
Compton Gamma Ray Observatory (CGRO) to be extremely luminous, hard 7-ray 
sourceJS^tSI -pj^g compact radio sources are members of a larger class of sources 
(constituting '--^ 10 % of all AGNs), designated as radio loud, that exhibit radio jets 
extending over many decades in radius (see Ref. 7 for a review) . The rapid variability 
and the superluminal motion of radio knots often seen in blazars indicate that those 
jets expand relativistically. According to the unified model both compact and 
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extended radio sources belong to the same class of physical objects, distinguished 
observationally by orientation to the observer; jets pointing in our direction are 
classified as blazars. 




Fig. 1. Schematic illustration of the high-energy emitting jet model. A relativistic jet is ejected by 
a central engine consisting of a black hole surrounded by a hot accretion disk. Soft X-ray photons 
emitted from the inner disk radii (indicated by the dashed lines) may be intercepted directly by 
the jet (red lines) and/or Thomson-scattered across the jet by matter surrounding it (black lines). 
External photons penetrating the jet interact with relativistic electrons and protons accelerated 
in situ, resulting in emission of high-energy 7-rays and very high-energy neutrinos. Synchrotron 
emission by the relativistic electrons accelerated in the jet contributes additional source of seed 
photons. As the jet propagates through the ambient (interstellar of intergalactic) medium it drives 
a strong shock that decelerates it. The shocked ambient gas forms a hot cocoon that gives rise to 
additional emission from much larger scales. The radio lobes seen in FR2 sources and microquasars, 
and the afterglow emission seen in long GRBs are clear diagnostics of the shocked ambient gas. 
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Nearly a hundred blazars (66 high-confidence identifications and 27 potential 
identifications) are listed in the 3rd EGRET catalog*^ as 0.1 - 20 GeV 7-ray sources. 
Their distances and 7-ray luminosities span a wide range, with the most powerful 
sources (e.g., 0528-1-134, 4C38.41) exhibiting isotropic equivalent 7-ray luminosities 
as high as 10"*^ erg s,~^ . Despite observational efforts no extended radio sources 
or radio-quiet AGNs have been detected at 7-ray energies (except for a marginal 
detection of Gen A) . This apparent exclusive association of EGRET AGN sources 
with compact radio sources strongly supports the unified model, suggesting that 
the 7-rays are produced inside the jet and are beamed, with a beaming factor 
/h ~ 0^/2 ^ 10^"^ — lO^'^, where is the jet semi-opening angle. The energy spectra 
in the range 0.1-10 GeV can be fitted by power laws, with energy spectral indices 

in the range 0.7 to 1.4, (though the broad energy response of the instrument 
coupled with the paucity of photons cannot really exclude spectral curvature) . 

Some blazars also exhibit very high energy (VHE) 7-ray emission. To date about 
a dozen blazars have been reported as VHE (> 100 GeV) source The energy 




Fig. 2. A sketch of a microquasar system. 
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spectra above 100 GeV are best fitted by power laws with upper cutoff at around 
a few TeV. Essentially all of the VHE sources are high-frequency-peaked BL Lac 
(HBL) objects located at relatively low redshifts {z < 0.2). The only exception is 
the FRI radio galaxy MSlEI. The BL Lac objects are blazars that show very weak 
emission lines or no emission lines at alQ, in contrast to the flat-spectrum radio 
quasars (FSRQ) that exhibit broad optical emission lines. This suggests that in BL 
Lac objects disk emission is strongly suppressed. 

The broad band uFi, spectrum of blazars is characterized by two main spectral 
components. A low energy component peaking in the submm-to-IR regime in FSRQs 
and in the UV-to-hard X-ray band in HBLs, and a high energy one peaking at MeV 
energies in FSRQs and GeV energies in HBLs. The ratio of luminosities of the high- 
and low-energy spectral components is typically large in the powerful sources, and 
is of order unity in the BL Lac objects (but may change during outbursts). The 
beamed radio-to-UV emission is most likely synchrotron radiation by non-thermal 
electrons (the synchrotron spectrum may extend up to hard X-ray energies in BL 
Lac objects). The origin of the high-energy (hard X-ray/7-ray) emission is still 
debatable. It may be due to inverse Compton scattering of some seed radiation by 
electrons accelerated in situ, or it may have hadronic origirffSl. Conceivable sources 
of scattered photons include the (beamed) synchrotron radiation produced in the 
jet (SSC mechanism; e.g., Refs. fTT]-[l4|) or external radiation (ERG mechanism; e.g, 
Refs. ll5H19p . presumably continuum radiation emanating from the inner parts of 
an accretion disk, that is either directly enters the jet or, alternatively, reprocessed 
and/or scattered across the jet by gas surrounding it, as illustrated in Fig.[TJ 

Strong, rapid variability over the entire electromagnetic spectrum is one of the 
characteristics of blazar emission. Large amplitude variations of the 0.1-10 GeV 
ray flux on timescales of days or weeks are typical to many EGRET blazars 
Dou bling t imes as short as a few hours have been reported for some of the strongest 
flareJM. 

It should be noted that the variability time that could be measured at 
GeV energies at the time was limited by the sensitivity of EGRET, so it is con- 
ceivable that the high-energy emission in blazars varies on even shorter timescales. 
Indeed, faster flux variations has been observed in the TeV blazars. In the case of 
Mrk 421, flux doubling time as short as 15 min has been reportecPSl. This rapid vari- 
ability of the hard 7-ray emission provides severe constraints on the bulk Lorentz 
factor of the outflow, and on the compactness and location of the emission region. 
It also has important consequences for the production of VHE neutrinos that will 
be discussed in f|6l 

Radio observations of EGRET AGN sources have shown that in man y cases 
radio outbursts seem to follow 7-ray flares with time lags of weeks to monthj2 5 | 26 | 
In a one or two exceptional cases, however, sources which have been monitored 



GeV 7- 
rj™. 



^Although M87 is believed to be a miss-aligned BL Lac the origin of the VHE 7-rays is unclear, 
as the emission from the jet is expected to be beamed away from us 

''Technically, a blazar is defined as BL Lac if the line equivalent width is less than 5 Angstroms 
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at radio frequencies for long time showed no radio activity following an EGRET 
flare. In general, time lags between 7-ray and radio outbursts are expected if the 
7-ray emission region is located well within the radio photosphere (see Fig. |6]) . 
Nonetheless, more complicated relationship can be envisioned. For example, in the 
radiative front model the delay between radio and 7-ray outbursts can be controlled 
by changing the opacity at the creation radius. In that case, one may expect to have 
some correlations between the ratio of amplitudes of the variations and the delay 




The relation between the 7-ray and optical/UV fluxes is not well characterized 

yet. In some sources the optica l, UV, X-and 7-ray fluxes seem to be correlated given 

the temporal resolutiorPSESl \yi the quasar 3C279 there are some indications for 

00 

lags of 1 to 2 days between the UV and 7-ray peak- , but these are not statistically 
significant. Future multi- wavelength observations may be able to resolve such short 
time lags, and provide additional constraints on the models. Observations of the BL 
Lac Mrk 421 reveal simultaneous X-ray/TeV outbursts with no significant variations 

on on 

of the flux at EGRET energies . Such an event has no natural explanation 
within the framework of inhomogeneous ERG rnodels. It may be reproduced by the 
one-zone SSC model under certain assumptioni^. 



2.2. Microquasars and 'y-Ray Binaries 

Microquasars are Galactic X-ray binary (XRB) systems, which exhibit relativistic 
radio jets 32-36^ These systems are believed to consist of a compact object, a 
neutron star or a black hole, and a giant star companion. Mass transfer from the 
giant star to the compact object through the formation of an accretion disk and 
the presence of the jets make them similar to small quasars, hence their name 
"microquasars." The analogy may not be only morphological; although there is 
no obvious scaling, it is common thinking that the physical processes that govern 
the formation of the accretion disk and the ejection of plasma into the jets are the 
same for both systems. Local, Galactic microquasars may therefore be considered as 
nearby "laboratories," where models of the distant, powerful quasars can be tested. 

There are, nonetheless, some important environmental differences that can affect 
the resulting high-energy emission from the system. In particular, in microquasars 
with a high-mass stellar companion (HMXBs) the radiation and the stellar wind 
emanating from the companion star can interact with jet, leading to additional 
emission of 7-rays and neutrinos. The compact object in some of those HMXB 
systems may be a pulsar. In this case high-energy emission may result from the 
interaction of the pulsar wind with the stellar radiation. Those systems are termed 
now 7-ray binarie^^. The emission from HMXBs is expected to be modulated, 
owing to the rotation of the binary system. 

The temporal behavior of microquasars appears to be rather complex. They 
exhibit large amplitude variations over a broad range of time scales and freq uencies, 
with apparent connections between the radio, IR, and soft/hard X-ray fluxes 
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The characteristics of the multi-waveband behavior depend on the state of the 
source, that is, whether the source is in a very high, soft/high or low/hard state^^. 
The ejection episodes are classified into several classes according to the brightness 
of synchrotron emission produced in the jet and the characteristic time scale of the 
event xhe duration of major ejection events (class A) is typically on the order of 
days, while that of less powerful flares (classes B and C) is correspondingly shorter 
(minutes to hours). The correlations between the X-ray and synchrotron emission 
clearly indicates a connection between the accretion process and the jet activity. 
Whether radio and IR outbursts represent actual ejection of blobs of plasma or, 
alternatively, formation of internal shocks in a quasi-steady jet is unclear at present 
(cf. Ref.|43|). 

7-ray emission from microquasars has been predicted shortly after their discov- 
ery I 44 | 45 l g^^j-jy attempts to detect microquasars with EGRET yielded only upper 
limits ESI. Tentative identifications of two EGRET sources with the HMXBs LSI+61 
30^33 and LS SOOoSSl-were reported later. Both systems have been detected recently 
by TeV observatories, confirming the EGRET identifications^^. Recent observations 
suggest that LSI +61 303 is a 7-ray binary. Possible association of some unidenti- 
fied, variable EGRET sources with Galactic black hole systems has been noted in 
Rei.m 



2.3. Gamma Ray Bursts 

The observed isotropic equivalent energy Eiso of GRBs span the range 10^" — 10^"* 
ergs (see e.g., Refs. [51] for recent reviews on GRBs). A large fraction of this energy 
is released during the prompt emission phase that, in long duration GRBs lasts 
for about several tens of seconds. This phase is followed by a longer phase during 
which afterglow emission at X-ray-to-radio wavelengths is typically observed. The 
spectrum of the prompt emission can be fitted by a b roken power law , with a 
peak energy that appears to be correlated with E'isJSSI (referred to as Amati rela- 
tion), albeit with a large scatter (but cf., Ref. [54]), and is about 1 MeV in the most 
luminous sources. In many sources the spectrum extends well above the peak, up 
to an energy of ~ 100 MeV in some cases, indicating that the emission cannot be 
produced in an adiabatically expanding e* fireball, as originally proposed^. Com- 
pactness arguments further imply that in sources exhibiting emission up to ~ 100 
MeV the Lorentz factor must exceed 100 or so. The large values of Eiso observed 
and the indications for achromatic breaks in the lightcurves of afterglow emissiorP^ 
are suggestive evidence for coUimation of GRB emitting outfiows ^^ ' ^'^ '. The distri- 
bution of opening angles inferred from the break of the afterglow lightcurves peaks 



iroi 

ai ^ O.W^. Beaming corrections seem to significantly reduce the scatter in the 
Amati relatiorP^', suggesting that this scatter is predominantly due to the spread 
in opening angles of the GRB population. It has been proposed that the correlation 
is entirely due to viewing angle effects and that the true energy of long GRBs may 
be standard with Ej '--^ 10 erg 
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The composition of GRB jets is yet an open issue. The standard view until 
recently was that the explosion energy is converted to kinetic energy of baryonic 
shells that collide and drive internal shocks at radii r ^ 10^^ — lO^fcm, behind 
which the prompt emission is produced. Production of 

UHECrJSHM! and VHE 

neutrino^- is then expected under certain conditions. Alternative scenarios have 
been proposed (e.g., Ref. [64]) . in which the MeV peak is limited by pair production. 
The recent detections of an early, shallow afterglow phase and sharply rising, large 
amplitude X-ray flares has introduced some difficulties to the standard fireball model 
and renewed the interest in baryon free fireballs (see §§ 16.41 for further discussion). 



3. Basic Jet Physics 

Let p', p' , e', and h' = (e' + p')/p'(? denote the proper rest mass density, total 
pressure, total energy density, and dimensionless specific enthalpy of the fluid, re- 
spectively (henceforth, prime refers to quantities measured in the comoving frame). 
The stress-energy tensor takes the form: 

T"^ = p'hu'^ul^ + p'g"!^ + ^{F'^'^F^ - ^g^f^F^), (1) 

47r Att 

where w" is the four-velocity, F^^, = d^Ay — di,A^ is the electromagnetic tensor, 
and 5"^ is the metric tensor. The dynamics of the MHD jet is governed by energy 
and momentum conservation: 

T% = 0, (2) 

mass conservation: 

(P'"");a = 0, (3) 

and Maxwell's equations. Under the assumption that the MHD flow is ideal (i.e., has 
infinite electric conductivity) Ohms law yields the additional constraint, u^F^p = 0. 
On scales of interest to the high-energy emission the gravitational force can be 
ignored. Using Eq. ^ the energy flux can then be expressed as 

T°'' = p'c^h'Tu^ + 5''=, (4) 

where F denotes the bulk Lorentz factor of the flow, = /T the corresponding 
3- velocity, and 5*^ the k component of the Poynting flux. The infinite conductivity 
condition implies that the proper electric field must vanish, viz., E' = E + v x 
B/c — 0. If in addition the magnetic field is toroidal, viz., B = B^(j), in which case 
B J- Vp, then the Poynting flux reduces to the simple form S = B'^Vp/Air. For a 
conical jet this is certainly a good approximation well beyond the light cylinder 
Rl = c/VL, where the ratio of toroidal and poloidal magnetic field components 
satisfies jS^/Bpl ~ R/Rl » 1, with R being the cylindrical radius. 

The total jet power is obtained by integrating T"*^ over a surface S perpendicular 
to the jet streamlines. Assuming a conical jet with an opening angle 9 and using 



5, 2008 14:27 WSPC/INSTRUCTION FILE jet-sub2 



High-Energy Aspects of Astrophysical Jets 9 

Eq. (g]) gives, 

= / T°'=dEfe = Tu''p'c^h'7r0\^{l + a). (5) 

Here a = B'^ / Anpc^h denotes the ratio of Poynting and kinetic energy fluxes, and 
is related to the comoving Alfven 4-speed through = c^/a (specificahy, is 
the component of the Alfven 4- velocity along the direction of wave propagation, as 
measured in the jet frame). Equation ^ implies that Lj is conserved. The mass 
flux carried by the jet is obtained by integrating Eq. 

M = p'u^'ne^r^. (6) 

This mass flux should be conserved to a good approximation, except for situations 
in which the rest mass content is significantly altered by processes involving, e.g., 
pair creation and annihilatiorEI or neutron leakage*^'''. Assuming M is conserved and 
setting a = and /i = 1 in Eqs. ([5]) and ([6]), it is readily seen that the asymptotic 
Lorentz factor of the jet is limited to. 

Too < (7) 

In the case of magnetically dominated jets for which a >> 1, acceleration to 
r = Lj/Mc^ requires effective conversion of magnetic into kinetic energ y. This 
may be difficult to accomplish through direct conversion in ideal MHD 
To illustrate this last point consider a magnetically dominated, conical jet sec- 
tion above the light cylinder, where the magnetic field is nearly toroidal. Magnetic 
flux conservation then implies v^B^ = u^B'^ cx r~^, which in turn implies that 
in the relativistic limit {v"^ ^ 1) the Poynting power Lb — T'^v^ [B'^ / ATi)TT9'^r'^ is 
to a good approximation conserved. Now, because Lj and Lb are conserved so do 
Lk = Lj — Lb and <j = LB/Lk, suggesting that at least on conical streamlines 
direct conversion of magnetic to kinetic energy is highly inefficient (see e.g., Ref [69l 
for a more detailed treatment). Similar behavior has been observed also in self- 
similar models. It has been argued that the slow conversion may be a consequence 
of the magnetic field geometry invoked. Indeed, recent analysis of the trans-field 
equation has demonstrated that effe ctive magnetic field conversion can be achieved 
in the case of a paraboloidal fielcP^. Other solutions to the sigma problem include 
magnetic reconnection and dissipation by MHD and plasma instabilities. 

Table 1 summarizes the range of absolute jet luminosities, Lorentz factors and 
variability timescales At inferred for different classes of high-energy sources. It is 
generally expected that certain properties will scale with the mass of the black hole, 
henceforth measured in units of solar mass, mBH = Mbh /Mq. In order to elucidate 
this scaling we find it convenient to measure luminosities in units of the Eddington 
luminosity, LEdd = IO^^ttibh erg s""'^, and distance in units of the Schwarzchild 



•^This includes pair cascade models (e.g., Refs. [191 1651 and neutrino-anti neutrino annihilation in 
GRB models (e.g., Refl66)l 
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radius, rg — 3 x IO^ttibh cm. We denote by £j = Lj/LEdd the dimensionless jet 
power and by f — r/vg the dimensionless radius. A canonical choice of parameters 
is £j ~ 1, rriBH ~ lO'* for AGNs; £j ~ 1, tubh ^ 3 for microquasars; £j ~ 10^^, 
TTT-BH ~ 3 for GRBs; and £j ~ 10^, uins ~ 1 for giant flare events in magnetars. In 
terms of these parameters the proper jet energy density is given by 

^'i = (ra\2 2 = 10'' erg cm'^, (8) 

and the proper magnetic field by 

B' . ^ 4 X 10« f^V%-(.r)- Gauss, (9) 

where ^b = Lb/Lj = a /{I + u). The comoving baryon density can be expressed in 
terms of 77p, the fraction of the jet energy carried by non-relativistic baryons, as 

< - 10'^ ""'-i',, (10) 



Table 1. Characteristic parameters of relativistic jet 
sources. 

Lj r At 

(org/s) 

GRB 10"'^-10=" 10^ - 10^ millisec - min. 

AGN 10"'2-10-*7 5 - 50 hours - years 

MQ 10^^-10**" 1 - 10 days 

OF lO^^'-lO**' 1 seconds 



3.1. Energy Dissipation 

Dissipation of the bulk energy occurs on various scales and by different mecha- 
nisms. The rapid variability of the 7-ray emission observed in blazars and GRBs 
suggests that a considerable fraction of the bulk energy dissipates already at rather 
small radii, 10^ — lO^r^. The dissipation mechanism may involve internal shocks, 
collimation shocks, and/or magnetic field instabilities. The latter possibility is par- 
ticularly relevant to magnetically dominated outflows in which dissipation by MHD 
shocks is inefficient. The conversion of magnetic to kinetic energy in Poynting flux 
dominated jets is a long standing issue. Models of shocks with ma gnet ic dissipation 
have been developed and used to calculate lightcurves in blazars.L^ However, the 
physics of magnetic field dissipation has not been addressed in those models. At- 
tempts to study particular mechanisms, including dissipation driven by magnetic 
reconnection, kink and toroidal field instabilities are presented in e.g., Refs. ,68^ l72l 
Quite generally, the radius at which waves created by a fluctuating source 
steepen into shocks is rdiss ^ r^7lc(5t, where 7^1 is the Lorentz factor associated 
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with disturbance speed (e.g., the Alfven speed), as measured in the fluid rest frame, 
and St is the characteristic timescale over which the parameters of the injected fluid 
varjESI. Since St ^ rg/c, we expect dissipation of the bulk energy to occur at radii 

fdiss > or, 

rrf,,, > 3 X lO^msz/r^ cm. (11) 

3.2. Blast-Waves, Cocoons and Afterglow Emission 

Additional dissipation occurs at radii lO'' — lO^'^^r^ behind so-called external shocks 
that result from collision of the jet with the ambient medium (see Fig. [T]). The 
radio lobes observed in radio galaxies (FR2 sources in particular) and microquasars, 
and the afterglow emission in GRBs are clear signatures of the associated blast 
waves. In GRBs and other type II supernovae the jet must first break out of the 
dense envelope of the progenitor star. This drives a radiation mediated shock that 
propagates outwards through the stellar envelope. Under certain conditions high- 
energy photons and VHE neutrinos may be produced during the shock breakout 
phase-4. 

The dynamics of the blast wave and the properties of the cocoon emission depend 
on the parameters of the ambient medium. In situations of impulsive energy injection 
as in, e.g., GRBs and other core collapse supernovae, the evolution of the system 
at sufficiently late times may become self-similar. In the ultra-relativistic case one 
can show that during the self-similar phase the Lorentz factor of an adiabatic blast 
wave declines with observer time as F oc t~I^J^ (see, e.g., Ref. [75| . This results in a 
power law decay of the late time afterglow flux in GRBs. Observations confirm that, 
at least in some cases, the self-similar model is a good description of the blast wave 
dynamics (see Ref. 1761 and references therein). However, the constraints imposed 
on the parameters of the shocked plasma are puzzling. Specifically, the strength of 
the post shock magnetic field inferred from observations is a factor ^ 10^ larger 
than the typical interstellar magnetic fielci^, suggesting that the post shock field is 
somehow generated in the shock and reaches MHD scales before decaying. Moreover, 
as mentioned in §§ l2.3l and discussed further in §§ 16. 41 recent observations of the early 
afterglow emission are inconsistent with the predictions of the self-similar afterglow 
model. 

3.3. Particle Acceleration 

The coUisionless shocks created in the jet, and the shocks driven by the jet into 
the surrounding medium, provide sites for particle acceleration. Both electrons and 
protons accelerate at the shock front. The mechanism commonly invoked to explain 
the inferred energy distribution of the non-thermal particles is first order Fermi 
acceleration This mechanism naturally produces a power law distribution of 
particles, dn/de oc e~^, with the exponent p depending on the shock compression 
ratio. For adiabatic shocks, test particle theory predicts p = 2 in non-relativistic 
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shocks and p ~ 2.2 in relativistic shocks in remarkable agreement with the 
observations. Energy losses and nonlinear effects may be important under some cir- 
cumstances and can alter the resultant spectrum. The diffusive shock acceleration 
theory does not address the so-called injection problem, namely what fraction of the 
thermal particles is injected into the Fermi process. The observations indicate that 
an appreciable fraction of the bulk energy is tapped for the acceleration of particles 
to non-thermal energies, suggesting efficient redistribution inside the shock transi- 
tion layer. The mechanism responsible for the scattering of particles upstream and 
downstream is yet another open issue. The scattering is most likely anomalous, 
involving collective plasma instabilities that are poorly understood at present. Fur- 
ther progress in our understanding of the microphysics involved requires detailed 
plasma simulations of coUisionless shocks. Nonetheless, some basic results that are 
directly relevant to the high-energy emission can be derived in a simple manner, as 
we now show. 

The maximum comoving energy of accelerated particles is determined by equat- 
ing he acceleration time with the relevant loss time. The time it takes to accelerate 
a particle to energy e' can be expressed as tacc = ^Wl/c, where = e'/e_B' is 
the corresponding Larmor radius, and 77 ^ 1 is an efficiency factor that depends on 
detailed physics. For protons radiative losses are typically small, and the relevant 
loss time is the escape time from the system, which in the jet rest frame is given 
roughly by t^sc — rO / cT . Equating tacc and tesc, using Eq. ^ and adopting ri — 1 
yields a maximum comoving proton energy, 

4.™a. < W/ = OeB'r ~ 3.6 x IQ^^isC^mBHY'^T-^ eV. (12) 

As shown in ii l5.3l below. at sufficiently small radii the maximum proton energy may 
be limited by photohadronic losses rather than escape and may be smaller than the 
confinement limit given by the last equation. 

The electron spectrum injected at the shock front is likely to be modified by 
radiative cooling, and conceivably pair cascades in sufficiently compact regions'^. 
The maximum electron energy is determined by equating the acceleration rate t~^^ 
with the total cooling rate. Including both synchrotron and inverse Compton cool- 
ing, and using again Eq. ([9]), the upper cutoff energy of the electron spectrum can 
be expressed as: 

e„c.l0^-5 1 + ^ (erf)V2 eV, (13) 

where denotes the comoving energy density of the soft radiation field inside the 
jet. 

The energy of thermal electrons depends on the shock velocity and the rate 
at which energy is transfered from ions to electrons. The mean energy per proton 
behind the shock is ~ (Fs — l )TOpC^, where Fg is the shock Lorentz factor as measured 
in the frame of the upstream fluid. The mean electron energy is some fraction of 
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the latter, which we parametrize a^ 

ee,th ~ (14) 

The thermal electrons just behind the shock will loose their energy quickly if the 
synchrotron cooling time, tsyn — Qirrnl^c' / GT^e.thB'^ , is much shorter than the 
dynamical time t^t = r/cF. Rapid cooling of thermal electrons (i.e., tgyn/td << 1) 
is expected at radii 

f<10-^. (15) 

Adopting Cj = 1, O^T^ ^1 — 10 for AGNs and microquasars, we conclude that 
fast cooling occurs at radii f ^ (10^'^ — 10^'^)£,e£,B ^ 10^ — 10^ in those sources. 
In GRBs, for which Cj = 10^^, fast cooling of electrons is expected on even larger 
scales. 

The evolution of the electron (and proton) spectrum is governed by a set of 
kinetic equations that incorporate particle injection and acceleration, particle cre- 
ation and annihilation, and relevant energy loss processes (see e.g., Refs. [311 [7^ . 
A complete treatment of the problem is rather complicated and requires detailed 
plasma simulations. Some parametrization of the microphysics is commonly invoked 
in order to simplify the analysis. The simplest approach is to assume that the elec- 
trons are, on the average, continuously injected and re-accelerated in a dissipation 
layer of width A ;$ r. If the cooling length is shorter than A a roughly steady elec- 
tron spectrum will be formed in the dissipation layer. Assuming impulsive injection 
and ignoring pair cascades, the electron energy distribution below the thermal peak 
(ee < £e,th) reads 

= K(r)7g , (16) 

where 7e = ee/weC^, and K(r) is some fraction Xe of the total electron density of the 
shocked plasma, that depends roughly on the ratio of the cooling length and A. If the 
jet composition is dominated by electron-proton plasma, then the electron density 
behind the shock can be estimated using Eq. ([5]), whereby we obtain (adopting 
(7 < 1 in the dissipation region) 

^(r) ~ ^£f^ = ioi9^^£f^(0Ff)-2 cm-3. (17) 

The factor Fs accounts for compression behind the shock. The energy distribu- 
tion above the peak depends on the injected spectrum. Assuming efficient injec- 
tion with a spectrum Qmj cx 7^*^ that joins smoothly the thermal peak, the elec- 
tron distribution above the peak can be written as dne(r, 7e)/d7e = ne{r)j^~j^^j~P; 

le,th < 7e < Imax, with p = g + 1. 



"^with this parametrization = 1 defines equipartition for Fs — 1 ~ 1 
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4. Radiation Processes 

4.1. Low Energy Emission and Target Radiation Fields 

4.1.1. Synchrotron emission 

The local synchrotron spectrum can be well represented by a broken power law, 
characterized by several break energies that depend on radius. At sufficiently low 
frequencies the synchrotron source becomes self-absorbed. The turnover frequency 
v'^^ at which the synchrotron absorption depth is unity, depends on the electron 
density at the corresponding energy. As shown below, the turnover frequency is typ- 
ically smaller than the peak frequency produced by the thermal electrons. Adopting 
for illustration the electron distribution given by Eqs. (TTBl) and pT| . which is ap- 
propriate in the fast cooling regime, one finds 

K„~3xlOM^^2^ (0r)-4/3r~-i eV. (18) 

\^ TTlBH J 

The comoving spectral energy distribution of synchrotron photons peaks at a fre- 
quency t-pg^j. = ll^th^'g, where -fe,th = (.e^thlra^c^, with ee,t/i given by Eq. HI]), 
is the random Lorentz factor of the thermal electrons, and 27rzyg = eB' /mf.c. By 
employing Eqs. ^ and (fT4|) we arrive at: 

/iz.;,,, 4 X 10^ (^-^\"\l[eTfr^ eV. (19) 

Note that the ratio ^'sm/^'peak °^ ie'^{^j''^BH/S,BY^^ is independent of radius and 
depends only weakly on the other parameters except ^g. It is therefore generally 
expected that z/^,,, << i^pg^^, for > 0.1. Note further that for GRBs, with niBH = 
3, ^bCj = 10", = 0.3, r = lO^Ea, 9 = 0.16'_i, and r = E^cAi adopted, where 
At = At_3 milliseconds being the observed variability time of the prompt emission, 
the observed peak energy is hupeak = ^hv'^ ~ 2 MeV {9-iT\/S.t-.'i)^^ . This is 




Fig. 3. The local spectral energy distribution of synchrotron photons, hv'n'^y^. The peak fre- 
quency i^'p^^f. and the turnover frequency u'^^ are indicated. 
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consistent with the observed peak of the prompt GRB spectrum, ahhough large 
scatter is expected for reasonable variations of the source parameters. 

The density of synchrotron photons per logarithmic frequency interval, 
n'gy„{r,v') = dn'gy„{r,u')/dhiu' , is shown schematically in Fig. [S] Above the 
turnover frequency n'^y„{r,v') oc (j^'/^peafc)""' with a = 1/2 below the peak and 
a — (p — l)/2 ^ 1 above the peak. In terms of the fraction £^syn — "syn/^ji where 
u'^yn — J hv'n'^y^dlniy' is the comoving energy density of synchrotron photons, we 
have 

<„(r, = 10^1 (t^) c--'- (20) 



where A = 2- 2{vUvpeakY'^ + 2[(i.;„, _ l]/(3 -p) is a normaliza- 

tion factor. It should be emphasized that for an inhomogeneous source the observed 
spectrum may differ from the intrinsic local spectrum. In particular, for an unre- 
solved conical jet the observed spectrum should appear flat, S^, oc (see, e.g., Ref. 

uni). 

The comoving density of synchrotron photons is constrained by the observations. 
For a source at a redshift z observed at viewing angle On (defined as the angle 
between the jet axis and the sight line), the observed peak frequency is Vpeak = 
<5i/pg^^/(l + z), where 5 = [r(l — /?cos6'„)]~-^ is the associated Doppler factor. The 
flux density at Vpeak, as measured on Earh, can be expressed as 

_ hcn'gynjrpeak) ( rpeakO \ r-k fr,^\ 

where is the luminosity distance, and rp^^ak is the radius from which photons at 
the peak originated, and is related to I'p^ak through Eq. p9|) . The index k depends 
on whether the photosphere is at rest with respect to the observer (continuous jet), 
in which case A; = 2, or whether it is moving (emitting blob) in which case fc = 3. 
Eq. (|21|) is of little use if the source is unresolved, which is usually the case at 
scales relevant to the high-energy emission. However, the rapid variability observed 
during intense flare states can be used as additional constraint on the source size. 
The latter is related to the observed variability time At via Vp^ak ~ cAtTS/{l + z). 
As seen, the comoving density of synchrotron photons inferred from the observed 
flux density is sensitive to the assumed Doppler factor; n'^y^ oc S^S~^'^^^\ 



4.1.2. Disk radiation 

An important source of ambient photons is the accretion disk surrounding the com- 
pact object. Denoting by Cd — Ld/LEdd the dimensionless luminosity of the accre- 
tion disk radiation, and assuming blackbody emission from the inner disk regions 
yields a rough estimate for the source temperature: Td ^ lO'^ {Cd/rnsHy^^ k. Con- 
sequently, the I'Fi, spectrum of the accretion disk radiation should peak in the soft 
X-ray band in stellar X-ray sources, and in the UV band in AGNs, in accord with 
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observations. However, the observations seem to imply the presence of additional 
components that contribute to the spectrum of the central continuum source. In 
particular, a power law extension of the peak emission up to hard X-ray/soft 7-ray 
energies is typically seen in both quasars and Galactic X-ray sources. This compo- 
nent is commonly attributed to a hot, tenuous corona above the disk that scatters 
disk photons. 

In Quasars the observed spectral energy distribution of the central continuum 
source is characterized by a big blue bump (BBB), with a peak energy that lies in 
the range hvBB ^ 20 — 50 eV. Above the BBB the spectrum can be approx imated as 
a power law with index s — 1.5 up to an energy of about 0.1-0.5 although 
it may slightly vary from source to source. At higher energies the spectrum flattens 
and has a slop of s ~ 0.5 up to the upper cutoff energy at about a few hundred 
keV. Below the BBB the spectrum flattens again to s < 1. A fraction £,d{r) of 
the disk radiation is intercepted by the jet at radius r. This includes both direct 
illumination and reprocessed emission, as shown schematically in Fig [T] (see ^j6] for 
further discussion). The associated target photon spectrum (photon density per 
logarithmic frequency interval) adopted below for blazars has the form 



n^^t{r, ~ 6 X 10 — — 

niBHr^ \25 eY J y^BB 



1/2 

1 



cm ^, 



(22) 

with s = 1.5 for v > ubb and s = 0.7 for v < vbb- 

The spectrum of the central continuum source is not as well constrained in mi- 
croquasars. Although extension of the thermal component to a few hundred keV is 
typically seen in those objects, it is not clear at present whether it is produced in 
the jet or in the disk corona. For our calculations below we shall adopt Eq. ([22|) 
with hvBB = 250 eV for the target photon density contributed by the external 
radiation source in microquasars. We emphasize that the pair-production and pho- 
topion opacities at most energies are determined by the portion of the target photon 
spectrum near and bellow the peak, particularly in microquasars. The seed photons 
above the peak interact mainly with the lower energy quanta. It is therefore antici- 
pated that the results presented below wont be sensitive to the shape of the target 
photon spectrum well above the peak. 



4.1.3. Stellar radiation field 

In microquasars associated with HMXBs, e.g., Cygnus X-1, LS 5039, LS I -1-61 303, 
the companion star, typically O or B spectral type, has a luminosity ~ 10^^ — 10'^^ 
ergs s~^, and a spectrum that can be well approximated by a black body spectrum 
with a temperature ~ 40000 k. The binary separation in those systems is of 
the order of a ~ 10^^ cm. At jet radii r ^ a the stellar photons incident into the 
jet at relatively large angles and are blue shifted in the jet frame. It is therefore 
naively expected that the external target photon field in the jet will be dominated 
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by stellar radiation at jet radii r > R{^dLdl L^Y^"^ ■, where = + (see Fig. [2]). 
However, the Thomson optical depth of the stellar wind may largely exceed unity 
at jet radii r ^ a (see Eq. [JS] below), in which case the jet may be shielded on 
these scales. Detailed analysis of Compton scattering of stellar radiation photons in 
microquasars is presented in Ref. [53] 



4.2. Inverse Compton Scattering: ERC and SSC emission 

The number density per unit frequency per unit radius of 7-rays produced by inverse 
Compton scattering of a soft photon distribution n'^ can be expressed as 



di^i^dr J "c^''^'^«'-7"-7 

where a'^{v'g, ^i'^) is a Compton scattering kernel that depends, in general, 

on the differential cross section and the energy distribution of relativistic electrons. 
In situations in which the seed radiation field in the jet is dominated by external 
photons (ERC emission) the comoving photon distribution is given by 

Ki^'s^tJ-'s^r) = {y'J^sfn^xt{vs.tJ-s,r), (24) 

where I/, = i/^r(l + /S^i^), /i. ^ (/i^ +/?)/(! + /^^f^), with /?= (1 - r-2)i/2 being the 
bulk 3-speed . For an isotropic power law distribution Uextii^sT IJ-s,f) — "■o('")'^7"; 
i^s,min < Vs < i^s,max, the comoving photon distribution is beamed in the backward 
direction and can be approximated as (e.g., Ref. ,19|) 

2a+2pQ+l 

(a + 1) 

This can be readily generalized to a broken power law or to the standard soft 
radiation field given in Eq. (j22[) . Expressions for the comoving soft-photon density 
in the case of direct illumination are presented in Ref. 84] The spectrum of scattered 
photons can be computed using Eqs. ((23)) and ([25)1 once the electron distribution 
is specified. Redistribution functions of inverse Compton scattering are presented 
in Ref. [19] for the standard soft spectrum adopted above. Equation (fT3| suggests 
that 10 GeV and 1 TeV photons can be produced in microquasars and blazars, 
respectively, even at the smallest radii (but they will be quickly absorbed if produced 
below the 7-sphere). From Eq. (|14|) we further conclude that in flat spectrum radio 
quasars and in microquasars thermal electrons scatter disk photons near the black 
body peak to observed energies hv-y ~ lO^F^^g/iz^BS ~ 0.1 — 10 MeV. This may 
be the origin of the MeV bumps seen in those objects. Further discussion on the 
applications of ERC models to individual sources is given in § [6] 

Synchrotron-self Compton (SSC) emission refers to situations whereby the syn- 
chrotron photons are Compton scattered by the relativistic electrons accelerated 
in sitrPJ^Tl^. In one-zone SSC models the same electron population producing the 
synchrotron photons scatter them to higher energies. Such models have been used 
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to infer the parameters of the emission zone in TeV blazars (e.g., Ref. I3ip . As men- 
tioned in § 12. li the basic assumption is that the synchrotron emission gives rise to 
the low energy component of the observed SED, while inverse Compton scattering 
of the synchrotron photons produce the high-energy component. The observed cut- 
off frequencies Vs and Vc of the low and high energy components, respectively, can 
be used to relate the magnetic field and maximum electron energy je,max'mec'^ to 
the Doppler factor 6. For j/g we have 

Photons having frequency Vc are produced by IC scattering in the Klein-Nishina 
regime when 

> ^ 1020 Hz. (27) 

This is typically the case in the TeV BL Lac objects. In that case we have 

Sje,maxmeC^ = hv^. (28) 

In sources for which the scattering is in the Thomson regime {vsv^y/^ « 10^° Hz, 
condition should be replaced by 

llmaxt^s = l^c- (29) 

For the TeV blazars Vs ^ Hz and i/c ~ 10^^ Hz during quiescent states. Eqs. 
(P7)) and ([25]) can be solved to yield 

7e.ma. = 10^-^<5-l(z.,/102^Hz), (30) 

5 = 5 X 10-^(5(iyc/10^^ Hz)"2(i/^/10^* Hz) Gauss. (31) 

In cases where the Doppler factor can be estimated the parameters of the emission 
region are fully determined (see §[J6?2]). 



4.3. Pair Creation and Annihilation 

At sufhciently small radii the compactness parameter, l-y = L^ot I {raeC^r) ~ 
10''£^f~^, largely exceeds unity for all classes of sources considered here. It is, 
therefore, expected that 7-rays will not be able to escape from the inner jet regions 
without creating pairs. Moreover, as will be shown below the bulk outflow energy 
cannot be carried by e^ pairs at small radii, owing to rapid pair annihilation. This 
implies that the energy released by the central engine must be transported to larger 
radii by baryons and/or magnetic fields. In what follows we consider 7-ray absorp- 
tion and pair annihilation in some detail. 
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4.3.1. The "f -sphere 

7-rays produced in the jet can combine with soft photons to produce electron- 
positron pairs. The pair production cross section is given by 



16 



(3-/34)ln(i±|^)-2/3(2-/?^^ 



(32) 



where /3 is the speed of the electron and the positron in the center of momentum 
framc^"^, and is related to the 7-ray energy e'^, the soft (target) photon energy hi/'^ 
and the cosine of the angle between their momenta, /i, through 



(l-/3') = 



2to^c'* 



0<^< 1; 



(33) 



(1 - ^J)hv'^e'^ ' 

(3 — Q, ^ — —\ corresponds to the threshold for pair production that we express as 



2.5 X 10 



11 (h^y 



VleV 



eV. 



(34) 



As a rough estimate one can adopt cr^-y ^ 0.2ctt near threshold. For a comoving 
soft photon distribution n'g{^,v' ,r) (per unit volume, per Ini^', per steradian) the 



7.6 - 




-3 -2 

log(e/l TeV) 



Fig. 4. Dimensionlcss 7-sphcric radius f-y = r-y/vg versus 7-ray energy e-y. Solid lines show f5^^* for 
^bbIvx = 0.1 and ^aCdihiy/25eV)-^ = (blue), 10"^ (red), lO"* (green) and 10"^ (black). 

The dashed lines show f^^" for {e,'^^synCj /Ae'^F^ 5) = IQ-^ (blue), 10"^ (red), 10""' (green) and 
10-^ (black). 
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e;, = 6 X lOM C'^rr eV. (36) 



pair production opacity at some radius r is 

J-l Jln[2/(l-p)e;] 

Both the external radiation intercepted by the jet and the synchrotron photons 
produced inside the jet provide an opacity to pair production. The corresponding 
optical depth can then be expressed as r^-y = r^i^* + r*^" where r^i^* denotes the 
contribution of external photons and t^^" the contribution of synchrotron photons. 

Consider first the contribution of the synchrotron radiation field. Denoting by 
^7 peak threshold 7-ray energy for pair creation by interaction with synchrotron 
peak photons, and using Eqs. ([19]) and p4|) one finds 

Assuming that the photon distribution is isotropic in the jet frame, and substituting 
n'g — n'gyn in Eq- (|35p. where n'^y^ is the distribution of synchrotron photons given 
by Eq. we obtain 

with a = 1/2 for e'^ > e'^^p^^k and a ^ [p - l)/2 ~ 1 for e'^ < e'^ ^^^^. For 
AGNs, GRBs and magnetars we anticipate r*^" >> 1 at e!!^ = ^'-y peak^ while for 
microquasars r^^" ~ 1. This implies that at any given radius only 7-rays having 
energies e'^ < ^^^^ can escape the system. From Eq. (1571) it is readily seen that the 
'7-sphere', defined as the radius r^^"(e-y) beyond which the pair production optical 
depth to infinity is unity, viz., r^^"(r^^", e^) = 1, is proportional to e'^. In terms of 
the observed 7-ray energy e-y — 5e'^ we have 

,,;.~ = 5xl0«(%^)(^), (38) 

where Eq. ([55)) and ([57|) have been used with a — 1. In the case of GRBs we conclude 
from Eq. ([55]) that 7-ray absorption by pair production may be important at radii 
f ;$ 10'' even at MeV energies (see e.g., Ref . for detailed calculations). 

To compute the optical depth to pair production on external photons, we trans- 
form Eq. (|35p to the star frame and take — next- We further assume that the 
external radiation field is roughly isotropic, which is a good approximation in situ- 
ations whereby the target radiation field is dominated by reprocessed or scattered 
disk photons. With these simplifications the corresponding optical depth reads 



T^^* ~ 2 X 10« 



^d^d\ ( hvBB 



\25 eY J V'=7,-BS 



1/2- 

1 



__]_ 



(39) 



where e-^^BB — 10GeV(/ii^i3s/25 eY)'^ and ,bb{vbbIvx)- The associated 

7-spheric radius can be readily obtained from Eq. ([55)) and the equation r^:^*(r, e^) = 
1. The 7-spheric radii rf,^" and r^^* are exhibited in Fig.|4| 



5, 2008 14:27 WSPC/INSTRUCTION FILE jet-sub2 



High-Energy Aspects of Astrophysical Jets 21 

4.3.2. Pair annihilation 

As argued in §§[33]in the inner jet regions the coohng time of electrons and positrons 
is much shorter than the dynamical time and, therefore, all pairs are expected to 
quickly cool to subrclativistic energy. The comoving pair density at a given radius, 
n'j. (r) , is then limited by electron-positron annihilation^. Equating the annihilation 
time, t'^nn ~ {o'ann'n'j.v±)~^ , wherc v± denotes the thermal velocity of pairs in the 
jet frame, with the outflow time t'j: ^ r/cT gives n± ~ Tc/aannV±f- The maximum 
power that can be carried by pairs at radius r is then 

L± ~ + n'_)T^m^c^TTe^r^ ~ 2r^meC^ne^r /aann- (40) 

Consequently, at radii smaller than the annihilation radius, 

faun ^ 10^^, (41) 

L± < Lj and so the jet power cannot be carried by pairs. For powerful blazars and 
microquasars this means that some alternative carrier of energy and momentum 
must be present at radii r ^ 10^ — 10^. For GRBs Eq. (^1]) gives farm ~ 10^. We 
note, however, that the above analysis does not hold in compact regions where 
the comoving pair temperature exceeds ~ 1 MeV, e.g., the base of a GRB fireball 
(see Ref. [55]l . In this case kinetic equilibrium is established whereby pair creation 
balances pair annihilation. However, as the pair outflow accelerates the comoving 
temperature quickly drops owing to adiabatic cooling and the pairs will eventually 
annihilate. This typically happens at radii f ^ 10^ — 10'^, well below fann- 



4.4. Inhomogeneous Pair Cascade Models 

Injection of electrons and positrons to energies well above the local threshold 7-ray 
energy would lead to a development of intense pair cascades. In the inhomogeneous 
pair cascade modeff^ the pairs are injected with high energy over many octaves 
of jet radius. At a given radius, 7-rays are produced through inverse Compton 
scattering of external disk photons. The energy of a freshly emitted 7-ray is then 
degraded, via pair cascades, into many 7-rays with energy near the threshold energy 
below which the 7-rays can escape to infinity without being further converted into 
pairs. This threshold energy increases with radius, as seen in Fig. [4l so that higher 
energy 7-rays tend to originate from larger radius. 

For sufficiently flat electron injection spectra, the resultant 7-ray spectrum re- 
flects the spatial profiles of energy dissipation and target photon intensity, because 
the asymptotic flux is dominated by 7-rays that emerged from near to their asso- 
ciated 7-spheres. For K{r) (X r~P in Eq. (fT6|l and £,d oc in Eq. ([22]) the emitted 
7-ray spectrum is roughly a power law with spectral index a-y ~ ap/(l -f g), where 
a is the average spectral index of the target photon fielcff^. 
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4.5. Constraints from Variability 

Constraints on the source parameters can be derived from the variabihty of the 
high-energy emission. Suppose that a variabihty timescale Ai has been measured 
at some observed 7-ray energy e-y. This imphcs that the emission at the observed 
energy originated from radii r ^ r^ar = T5clS.t/ (1 + z), where 5 is the Doppler factor 
and z is the redshift of the source. At the emission zone the pair production optical 
depth at the observed energy must not exceed unity, viz., T^^^{r^an^'j/ 5) < 1, 
implying r^^"(e^) < r^ar- Suppose now that the synchrotron flux near the peak 
can be measured simultaneously. Using Eq. ([2T|) with fc = 3 and the fact that 

fpeak < Tyar WC dcduCC that 

For a low redshift source [z < 1) we take ~ 10^^h~^z cm. By combining Eqs. 
(|20]) . ((38)) and (gS]), and adopting ^ = 5 we arrive at 

iee^syn/^By^'S' ^ 2 X lO'' At-'^^Ve)-' [e^ / ITeY)'/^ S , (43) 

where Sjy is the flux at i^peak in Janskys and Ath is At measured in hours. 

5. Hadronic Processes 

Relativistic jets are also potential sources of ultra-high energy cosmic rays (UHECR) 
and VHE neutrinos. This of course requires that a considerable fraction of the bulk 
energy will be carried by baryons. In contrast to electromagnetic emission that can 
have either leptonic or hadronic origin, VHE neutrino emission is a unique diagnostic 
of a hadronic content. High-energy neutrinos can be produced in astrophysical jets 
mainly through the decay of charged pious: 



+ e'^ + Ve + i^,. + (44) 



Decay of neutral pions. 



tt" ^ 7 + 7, (45) 

leads to production of high energy photons, and under certain conditions may com- 
pete with inverse Compton scattering. The pions may be produced through collisions 
of protons with target photons (p7jj, or via inelastic (pp) and (pn) collisions. The 
latter mechanism produces both 7r+ and tt~ in roughly equal numbers, whereas 
the former mechanism produces mainly 7r+ in regions of moderate optical depth. 
Therefore, it may be possible to discriminate between the two production modes of 
charged pio ns by m easuring flavor composition (e.g., Refs. [89H92|) . or even to probe 
new physicj^ ^ l I, although the interpretation of such measurements may not be 

'Ho' 

straightforward^ . 

"^In cases whereby the neutron lifetime is longer than the loss time due to photomeson interactions 
the production rate of via (ri'y) collisions becomes comparable to the 7r+ production rate. 
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5.1. UHECR Production 

As explained in § 13.31 a substantial fraction of the energy dissipated behind shocks 
can, in principle, be taped for acceleration of protons to ultrahigh energies with a 
power law spectrum dn'^/dtp oc e'j^^'^ and an upper cutoff t'p^rnax that is limited by 
confinement (see Eq. [T^ . Consequently, production of UHECRs of observed energy 
tp = Te'p requires 

(.BniBHC, > 10«(ep/1020-5eV)2, (46) 

which essentially leaves GRBs with ttibh ^ 3, ~ 0.1, Cj ~ 10^^^ and powerful 
blazars with itibh ^ 10^, £,b — 0.1, Cj ~ 1 as the main candidates for astrophysi- 
cal UHECRs source^. Whether protons can actually be accelerated by relativistic 
shocks to the highest energies observed is yet another issue. 

5.2. Inelastic Nuclear Collisions 

Pion production via interactions of the ultra-high energy nuclei accelerated in the 
jet and the cold jet material is typically inefficient in sub-Eddington sources, but 
may be relevant in highly super-Eddington sources. For the baryon density given in 
Eq. pUj) the optical depth for inelastic nuclear collisions is approximately 

r,,^a,Xf -10-^1^, (47) 

where app — 50 mbn has been adopted, which is suitable for rough astrophysical 
estimates. For Eddington sources like AGNs and microquasars, for which Cj < 1, 
this is always much less than unity. In GRB models that invoke r/p ~ 1 at F ~ 100, 
9T ^ 10 and Cj ~ 10^^ this process may be relevant at radii f ^ 10''. 

5.2.1. Target ions from stellar winds in microquasars 

In microquasars with high-mass stellar companio ns a pot entially important source 
of target protons is provided by the stellar 

wincPM, as illustrated in Fig. [2] 
Typical mass loss rates and terminal velocities of O stars are of the order of Af^ ^ 
lO~^M0yr^^ and Vw ~ 2000 km s^^, respectiveljJ^. At a distance R from the 
companion star the wind density is Up ~ / {ATTrripR^ v^) ^ and the optical depth 
for pp collisions is 

If stellar wind ions can penetrate the jet then effective production of pions in the 
jet/wind interaction site, and consequent emission of VHE neutrinos and 7-rays is 

^ Condition II46I I does not hold in systems that violate the ideal MHD condition, e.g., vacuum gaps 
in dormant AGNs^? §8. ^nd magnetars. In this case the luminosity is limited only by the rate at 
which UHECRs are produced. 
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anticipated. Such an origin has been proposed for the TeV emission from the two 
Tev microquasars detected thus far. Kelvin-Helmhohz instabihties at the interface 
between the jet and slow wind may give rise to entrainment of the wind material. 
Whether significant fraction of the stellar wind ions can be intercepted by the jet 
is yet an open issue. Owing to rotation of the binary system modulation of the 
observed fluxes is expected. 

5.2.2. Neutron pick up in GRBs 

In some scenarios, the ultra-relativistic GRB-producing jet is envisaged to be en- 
sheathed by a slow, baryon rich wind emanating from the disk surrounding the 
black hole. This wind contains free neutrons out to a radius of ~ 10^ — 10^^ cm, 
that can diffuse across magnetic field lines into the central baryon free jet. The leak- 
ing neutrons are then pic ked-up and converted to protons in a collision avalanchJ^. 
It has been suggested^ that baryon loading of GRB fireballs is accomplished by 
this process; the number of captured neutrons has been found to be in reasonable 
agreement with existing limits on the GRB baryonic component. The charged decay 
and collision products of the neutrons become ultra-relativistic immediately, and a 
VHE neutrino burst is produced with an efficiency that can exceed 0.5. Other signa- 
tures may include lithium, beryllium and/or boron lines in the supernova remnants 
associated with GRB's and high polarization of the 7-rayj22l_ 

5.3. Photomeson Production 

Collision of a proton with a photon can lead to neutrino and 7-ray production 
through the following reactions: 

p -I- 7 ^ n, -(- 7r+, (49) 
p -I- 7 ^ p + vr" . 

The photopion production cross section peaks at CTp-y ^0.5 mb at the A resonance 
and drops to ap~f ^ 0.1 mb at higher energies where multipion production dominates. 
The inelasticity is K^: ^ 0.2 near the threshold and K-j^ ~ 0.5 in the multipion 
production regime. The threshold proton energy for which head-on collision with a 
target photon is at the A resonance is 

4,*/.-7xl0i6(^^^ ' eV, (50) 

where hu'^ is the energy of the target photon, as measured in the jet rest frame. 

Consider now photomeson interactions on the target synchrotron photons pro- 
duced in the jet. For illustration we adopt the synchrotron spectrum given in Eq. 
pop . The comoving proton energy for which interaction with photons near the peak 
is at the A resonance is related to the corresponding pair production threshold 
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through Eqs. dMl), (IMl) and dSOl): 

^p,peak = 10'-'<,pea;c = 2 X lO^" ( ^] ' ^'HYf eV. (51) 

^ ^ \mBH J 

The energy loss rate of protons due to photomeson interactions is t~}^ ^ 
ap-fKT^n'gy^c. Equating the latter with the acceleration rate t~}^ ~ ei3'c/ep, us- 
ing Eq. (|20p with a = 1 and Eq. (jSip . and assuming single pion inelasticity factor 
of -ftTjr = 0.2 one obtains an upper limit on the proton energy: 



< 2 X 10^5 f-AV^V-^l '^'c'/'^R^'(^rr)i/3 eV. (52) 

\isynj \mBH J 



Comparing the latter limit with Eq. (jl2p then implies that the maximum proton 
energy is limited by photomeson losses at radii 



f <6xl03('Mf^Vf^') 



(53) 



and by confinement at larger radii. From Eqs. (fT2|) . ([51]) and (|52|) it is evident 
that protons can be accelerated to energies exceeding that required for photomeson 




Fig. 5. Dimensionlcss 7r-sphcric radius = rn/fg versus proton energy measured by an on-axis 
observer ep = re'p. The left panel corresponds to a choice of parameters that represents a typical 
blazar like 3c279: 9T = 1, T = 10, mgfj = 10*, and the right panel to a choice of parameters 
representing a typical microquasar like GRS 1915: S = 0.1, F = 3, rriBH = 3. The solid lines in both 
windows show fj^* for ^^Cdihi'/25eV)~^ = 10~^ (orange), 10"'^ (cyan), and lO""* (green). The 
dashed lines show rj^" for Cj = 1, = 0.1, = 0.3 and i,3yn/A = 10~^ (black), 10~^ (red) and 
10 (blue). The dotted-dashed lines correspond to the maximum proton energy ep^rnax = 
computed using the minimum of Eqs. I I12I I and 11521 . with the same color code employed for the 
dashed curves. The break in the dotted-dashed lines corresponds to the transition from radii at 
which the maximum proton energy is limited by inelastic collisions to radii at which it is limited 
by escape (see Eq. 1531 ). 
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interactions with synchrotron peak photons, viz., pg^j, < Cp max^ radii 

A rough estimate of the photopion optical depth contributed by synchrotron photons 
gives 



-"(4,r):.a,,n' J^O.lsfl^U^^) ( 1 , (55) 



1/2 



p^peak 



with a = 1/2 for > 6^^^^^, and a ~ 1 for < e'p p^ak- ^'^'^^ ^^'^ photopion 
opacity increases with increasing proton energy and is largest at e'p ^ax- This trend is 
expected for any reasonable synchrotron spectrum, at least up to a proton energy at 
which interactions at the A resonance involves synchrotron photons having energies 
above the self-absorption frequency hi^sm- analogy to the 7-sphere we introduce 
the TT-sphere, the radius r^^" that solves the equation t^'^ {e'p, rp^^'-) = 1. In terms 
of the transformed energy Cj, = Te'p we obtain from Eq. (j55p 

where a = I for 0.15(6ynM6'r2)(/:jmBH/Cs)^''^ > 1 and a = 1/2 when the 
opposite inequality holds. 

As for pair production, external photons intercepted by the jet also provide tar- 
gets for p7 collisions. The corresponding optical depth, computed using the external 
radiation spectrum given in Eq. (j22p . reads: 



P7 — 



\25eYj \ep,BB 



1/2-s 

I r ^ \ 

1 



(57) 



where ep^BB ~ 10^^'^eV(/ii'BB/25 eV)^^ and ep,x = ep,BB{vBB /i^x)- The 7r-spheric 
radius r^^*(ep) associated with the external radiation field can be readily obtained 
by setting Tp^* = 1 in the last equation. The radii f^^" and f^^* are plotted in figure 
[SJ A comparison of the different scales derived above is presented in Fig. [S] 



5.4. Radiative Cooling of Pions and Muons 

Radiative cooling of charged pions and muons can limit the energy of the resultant 
neutrinos and, more importantly, the overall efficiency of neutrino production. For 
a pion having energy — r7i^c^7^, the ratio of synchrotron cooling time and decay 
time, — 2 X 10^*7^ s, is given by 

^ = 5 X 10237;-2b'-2 = 3 X 106^;-2 f \9Trf, (58) 

\niBH J 

where B is given in Gaussian units. From the last equation it is evident that high 
energy neutrino emission is severely suppressed in regions where B > 10^^ G, such 
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as expected near the surface of pulsars and magnetars. The only exception is the 
gaps formed in charge starved regions in the magnetosphere, where the pions can 
accelerate along magnetic field lines by a parallel electric In general pions 

having energies 

e' > 2.4 X 10" ( ' QVf eV (59) 
\mBH ) 

will cool radiatively before decaying. For muons this limit is smaller by a factor 
[th/tt^Y/'^ = 10, where is the muon lifetime. 

Suppose now that protons can be accelerated to the maximum limit allowed by 
confinement, and consider the implications for synchrotron cooling of charged pions 
produced by photomeson interactions of those maximum energy protons. To eluci- 
date the limit on the maximum proton energy above which the resultant photopions 
lose their energy radiatively, we express equation ()58|) in terms of the confinement 
energy given in Eq. ((T^). Using Eq. ([55)1 and adopting inelasticity of = 0.2 for 
single photopion production; that is, taking — 0.2econf yields, 

tsyn 



2 X 10" 



(eco„//1020 eV)4"- ^^^^ 
This suggests that if protons are accelerated to energies in excess of 2 x 
lQ^^{OmBHfY^^ eV, then the consequent charged pions will lose most of their en- 
ergy before decaying. 

Pions also lose energy as a result of inverse Compton scattering of synchrotron 
photons. At pion energies 



> e.,c = 5 X 10^ 



msH 



-1/2 



evr 



eV 



(61) 



the scattering of peak synchrotron photons is in the Klein-Nishina regime. Com- 
paring this with Eq. ([?T|) we find that e^^^/^p peak — 0-25^^, which for '^-^ 1 equals 



f 



■^g ^diss 



r^'=P,max ) r^lTeV) r(lGHz) 




r(cm) 



Fig. 6. Characteristic scales exhibited for typical Galactic and extragalactic sources. The radius 
above which shocks can form r^i^g, the 7r-sphere at maximum proton energy ^'tt (£p,max), the 'y- 
sphere of a TeV photon r.^(l TeV), and the radio photosphere at 1 Ghz r(lGhz) are indicated. 
The shaded area corresponds to the region where effective production of VHE neutrinos can occur. 
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roughly the inelasticity factor ■ Consequently, for plasma at rough equipartition 
pions produced by protons of energy e'p > ^pp^ak satisfy > Cjr.c- The energy loss 
time due to inverse-Compton scattering is then approximately 



hv' 4^2 1 

. " -TT peak . t r\ — 5 >e^B / 

^--^ 8c..<,„ Kc2)2 -4x^» (62) 
where 7^ = e^/rn^c^, and the ratio of the latter and the pion lifetime is 

^^2x10^^^, (63) 

independent of pion energy. Since A ^ a few, it is concluded that inverse- 
Compton losses of pions and muons can be ignored provided ^l^s/Csyn > 
and (.lisl^syn > 10"^, respectively. 



5.5. Relations Between Pair- Production and 
Pion-Photoproduction Opacities 

Observations of VHE 7-rays can constrain the pion photoproduction opacity, par- 
ticularly in situations where rapid variability of the VHE 7-ray flux is observed. 
To illustrate the relationship between the pair production and pion photoproduc- 
tion opacities, consider first interactions of 7-rays and protons with the synchrotron 
spectrum presented in Eq. ((20|) . From Eqs. ((37|l and ([55| one finds 



T^r(e' r) 




Fig. 7. Opacity ratio Tp^" (e'p) / t^^" {e'^) as a function of comoving proton energy tp' at f = 10^ 
(dashed lines) and 10^ (solid lines). The blue, red, green and black colors correspond to comoving 
7-ray energies of e' = 10, 10^, 10"^, and 10'* GeV respectively. 
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where a = 1/2 if both > e^^ peafc ^'p > ^p.peak a = 1 when the opposite 
inequahties hold. A plot of this ratio in the general case is displayed in Fig. [T] The 
ratio of opacities associated with the external radiation field, Tpi^*(ep)/T^i^*(e^), can 
be computed in a similar manner and is displayed in Fig. [51 As seen from the figures, 
at 7-ray energies above a few TeV the opacity ratio is smaller than unity even at 
the maximum proton energy. 7-ray observations at these energies can therefore 
constrain the neutrino flux that can be emitted from the same region. A particular 
application of this result is discussed in §§ 16.21 



5.6. Neutrino Flux and Spectrum 

The neutrino flux emitted from the jet depends on several factors: (i) the fraction 
of jet energy injected as a power law distribution of protons ^p, (ii) the spectrum of 
accelerated protons, (iii) the cooling times of pions and muons and (iv) the pion pro- 
duction opacity. It is customary to introduce the parameter /^(ep) < 1, representing 
the fraction of proton energy converted to pions. In case of pion photoproduction 
it can be approximated by /^(ep) = niin[l, i4r7rTp-y(ep)]. For a proton distribution 



1 - 




log(E/lTeV) 

Fig. 8. Opacity ratio TpHj^ (ep) / t^^^ (ej) as a function of observer frame proton energy ep. The 
blue, red, green and black colors correspond to 7-ray energies of tj = 10, 10^, 10'^, and 10* GeV 
respectively. 
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'np{e'p) the average proton energy lost to pions can be defined as 

J ^ //7r(£p)£>p(ep)c?£; ^^^^ 

This average fraction depends on the spectrum of both protons and target photons. 
As argued in §§ 15.41 pions are expected to decay prior to significant energy loss, 
while muons may lose significant fraction of their energy before decaying. We shall 
therefore conservatively assume that a single high energy (or D^) is produced in 
a single photopion interaction of a proton (or neutron) , corresponding to conversion 
of 1/8 of the energy lost to pion production to muon neutrinos. In terms of and 
/tt the flux at Earth of and z/^ can be expressed as, 

^ ^J^r-'S'^ lO-^o^jJla^ erg s-hm-\ (66) 

L L'2S 

where S is the Doppler factor, and — W^^dL2& cm is the luminosity distance to 
the source. 

A rough estimate of the neutrino spectrum expected to be produced in sources 
of high optical depth to photopion production can be made using Fig [51 Firstly, 
assuming that neutrinos are produced at radii larger than the dissipation radius 
fdiss given by Eq. (jlip and recalling that the energy of the resultant neutrino is 
~ 5% of the proton energy, we expect significant v flux in the energy interval 30 - 10^ 
TeV in AGNs, 1 - 10^ TeV in microquasars, and above 10^ TeV for GRBs internal 
shockd^. Secondly, for a thick target the neutrino spectrum should be similar to 
the proton spectrum. Thus, we expect dn^/dei, oc gjj^ in the above energy intervals. 

The probability t hat a muon neutrino will produce a high-energy muon in a 
terrestrial detector is^EIlp^^ ~ 1.3 x lO-^E^.^^^, with P = 2 for < 1 TeV and 
P — 0.8 for El, > 1 TeV. Thus, for a flat neutrino spectrum above 1 TeV, the muon 
flux at the detector is « {Po/Eq)T„^, where Pq/Eq = 1.3 x lO^^TeV"^ The rate 
of neutrino induced muon detections is 



TV, ^ 2 X lO-^%fJ-l!!^M^iA/lkTn') day-\ (67) 

^ "L28 



where Eq. (|66)) has been used. Estimates of neutrino yields for each class of objects 
is given in § [HI 



6. Applications 

6.1. Flat Spectrum Radio Quasars 

The prodigious emission exhibited by the most powerful 7-ray blazars imply Cj > 1 
for a black hole mass range itibh ~ 10* — 10^" inferred for quasars. Lorentz factors 
F ~ 10 — 30, as inferred from superluminal motions, are typical for FSRQs. Line 
luminosities and observations of radio quiet (jetless) quasars suggest disk luminosity 
of Cd ^ 1. The main uncertainty is the fraction intercepted by the jet. As seen 
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from Fig. REGRET (0.1-20 GeV) 7-spheres are typically located within the broad 
emission line region. On such scales an important fraction of disk radiation can 
be Thomson scattered across the jet by the broad line gas. The fraction is then 
roughly equals the average Thomson depth < > of the scattering mediurrilS], xhe 
standard emission line clouds are generally modeled to have free electron column 
densities of ~ 10^^ cm^^ and covering factors ^ 0.1, giving tt >~ 10""'. However, 
the inter-cloud medium, which is often invoked as an agent to confine these clouds, 
is widely supposed to have an electron depth a factor ^ 100 larger (e.g., Ref. I104( ). 
Furthermore, bound electrons, which may have a much greater density than free 
electrons, will present a roughly Thomson scattering cross section to hard X-rays. It 
is therefore reasonable to expect 10"^- 10"^. Since (^g^syn-Cj/A^^F^J) < 10"^ 
is expected for the inferred values of F and Cj , it is concluded that the target photon 
field should be dominated by the reprocessed disk emission if ^ 10~^ (see Fig. 
HI) . Detailed models that take into account SSC and ERG component show that 
the ratio of ERG and SSG luminosities required by spectral fits lies in the range 1 
to 10. 

The 7-spheres can be mapped in principle by measuring temporal variations 
of the 7-ray flux in different energy bands during a flare. If the 7-ray emission is 
indeed produced over many octaves of jet radius where intense pair cascades at 
the observed energies are important, as in the inhomogeneous pair cascade model 
discussed in §§ 14. 4( then we expect the high energy to either lag the lower energy or 
to vary more slowly. With the limited sensitivity and energy band of the EGRET 
instrument it was practically impossible to resolve such effects. It is hoped that with 
the upcoming GLAST instrument this will be feasible 

As seen from Fig.O above the dissipation radius rdiss — F^ ~ 10^, the photopion 
optical depth exceeds unity at proton energies above 10^^ eV for £,dCd — 10~^. The 
maximum proton energy there exceeds lO'^^ eV. Assuming np(e'p) oc e'~'^ we estimate 

~ 0.5. Likewise, for ^dCd = 10"^ wc find ~ 0.2. With tubh = 10^ Cj = 1, 
F = 10, = 0.1 and ~ 0.2 Eq. ([57]) implies detection rate of about one signal 
event per year in a cubic-scale neutrino detector. A similar conclusion has been 
drawn in Ref. 11031 who specifically considered the flare observed in 1996 from the 
quasar 3G279. 

6.2. TeV Blazars 

As mentioned in 12. H all known TeV blazars are BL Lac objects, which are typically 
much fainter than FSRQs. The observed bolometric luminosity of the TeV blazars 
during quiescent states is of the order of a few times 10'*'* ergs s~*, with about 10 
percents emitted as VHE 7-rays. The VHE flux increases by more than a factor 
of 10 during flaring states. In terms of the observed luminosity of VHE photons, 
denoted henceforth by Lvhe, and the corresponding radiative efficiency ^rad-, the 
jet power can be expressed as 

= Cadi^Vry'LvHE. (68) 
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Estimates of the black hole mass in a sample of sever al dozens BL Lac objects 
yield a range of 10^ < ttibh < 10^'^, similar to quasar Jl^. In particular, for the TeV 
blazars Mrk 421, Mrk 501 and lES 2344+514, Ref.[105]estimates log{mBH) = 8.22± 
0.06, 8.62±0.11 and 8.74±0.18, respectively. This imphes Cvhe — Lvhe/ LEdd ^ 
10-3 _ 10-2. 

The absence of emission lines suggests, as often argued, that the target radiation 
field in the BL Lac objects is dominated by the synchrotron photons produced in 
the jet. At radio to IR frequencies the observed flux lies in the range 10-^ — 1 Jy. 
Adopting r6' = 1, Sjy = 0.1, and z = 0.03 (roughly the redshift measured for Mrk 
421 and Mrk 501) Eq. (gS]) yields 

5>lQ"\^B/ee^synf"°/^tl"''. (69) 



As described in § l2.1l above. variations of the VHE 7-ray flux on timescales At/j ^ 1 
have been reported for Mrk 421 and some other BL Lac objects. Since quite generally 
(Se^sj/n/Cs) ^ 1 Eq. ([69]) imphes 6 ~ 10—50 during flaring states (see also Refs. 11061 - 
llOSp . Such high value s are in c lear disagreement with the much lower valu es inferred 
from source 

statisticJininni 

and superluminal motions on parsec scaleJIUl^II^l, 
Various explanations, inclu ding jet deceleratior J-*^ ^ 5 | 1 1 6j^ a str ucture consisting of 
interacting spine and sheatlPSi, and opening ang le effectjnil 

have been proposed 

in order to resolve this discrepancy. Under the one-zone SSC model Eqs. (|5D|) . PT|) 
and (|69p give a maximum electron energy meC^Jcmax ~ 0.1 TeV, magnetic field 
B ~ 0.1 G, and Doppler factor (5 ^ 30 in the emission region at a radius of r ^ 10^^ 
cm. 

Consider now the implications for 7-ray and neutrino emission. First, with the 
above estimates of the black hole masses we have cAt/rg ^ lAth. Since the 7- 
spheric radius must be smaller than the characteristic size of the emission zone 
we deduce that f.y{e^ = ITeV) ^ 5"^. At small viewing angles 5 ~ 2r, implying 
f^(e^ = ITeV) ~ rdiss, where rdiss is the dissipation radius given by Eq. (jlip . 
Thus, the observed variability seems to imply that the unsteadiness of the flow 
parameters occurs over a dynamical timescale. Assuming the source is observed 
at small viewing angles, viz., 5 ~ 2r, and adopting Cvhe ~ IQ^^ (roughly the 
value measured for Mrk 421 and Mrk 501) in Eq. ^ gives C-j = lO-^CadE"^- 
Substituting the latter result into Eq. ifTS]) and taking 9T = 1, we obtain 
^e,max ~ 10^^'^ (^raci/Cs ) ^^"^ (E/IO)'^/^ eV, consistcnt with the observed spectrum. 
Likewise, Eq. ^ gives e^^rnax < W^^-^{£.B/£.radV^^{^/W)''^ eV. Combining the 
latter result with Eq. ^ we find ,„^^/e^ ^^^^ < 0.03^^ (^/10)-^ which combin- 
ing with Eq. (IMl) imphes Ep.mas/ep peafc < 1- ^^nce T^^"(e7 = ITeV) < 1 we deduce 
from Eq. (gl that r^^^" < 4 X IQ-^Cp/ep^^ax) and U{e'p) = K^t^v^' ~ 0.2rp^J;" < 
^Q~^{e'p/ f'p^max)- Substituting /^r into Eq. ((65|) and assuming np{e'p) cx e^-^ have 

< 10-3/ln(e^ „,„^/e^^„,„) < lO"*, where e^^„„ ~ IGeV is the lower cutoff en- 
ergy of the proton energy distribution. This demonstrates that, contrary to earlier 
claims, the emitted flux of high energy neutrinos should be only a small fraction of 
the observed VHE 7-ray flux. With the above numbers the rate of neutrino induced 
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muons expected to be measured by the upcoming detectors is 

< 0.03ep(A/lW) yr-\ (70) 

too small to be detected. 
6.3. Microquasars 

In microquasars with resolved jets a minimum jet power can be estimated from 
the observed brightness temperature of radio kn ots. Such estimates yield Lj ~ 
10'^^ — 10'^^ ergs s^^ for most source j^ ^ l 1 1 8 1 1 1 tH p^j. gyg^g^is with a measured 
black hole mass the jet power appears to be close to the Eddington limit, viz., 
Cj ~ 0.01 — 1. The contributions of the jet, disk and companion star to the ob- 
served broad-band continuum spectrum cannot be directly extracted from the obser- 
vations. It is believed that in microquasars associated with low-mass X-ray binaries 
the accretion disk is the dominant source of external photons, whereas in micro- 
quasars with a massive stellar companion t he stellar ra diation makes an important 
contribution to the external radiation fielc|S2llllIl2ni, xhe observed X-ray emis- 
sion, whether due to disk emission or synchrotron emission from the jet implies 
r-f » Tdiss at GeV energies, suggesting that, like in blazars, pair cascades may be 
important. Detailed treatment of 7-ray emission from leptonic microquasar jets and 
fits to individual sources are given in e.g., Refs. l44l [45l [83l 1120111221 Hadronic 7-ray 
emission has also been considered (e.g., Ref . l94l [95]) . This is particularly relevant 
for microquasars with high-mass stellar companions for which, as can be seen from 
Ea. (|48)) . neutral pions can be effectively produced by inelastic colhsions of relativis- 
tic protons accelerated in the jet and the ions in the stellar wind. In this case the flux 
of VHE neutrinos should be comparable to the observed 7-ray flux. Neutrino yields 
have been estimated for the TeV microquasar LS I +61 303 assuming the observed 
TeV emission has a hadronic origirPSl. With the upcoming neutrino experiments it 
would be possible to constrain this model. 

Next, we estimate detection rates for VHE neutrinos produced by decay of 
photopions in intermittent microquasars (see Ref. 11231 for more details). The 
latter exhibit ejection episodes that lasts several days, with a total energy of 
Ej = LjAt ^ 10^^ — lO''^ ergs. The expected number of events in a single burst is 

7V^^20Cp/.^^(^/W), (71) 

where Ej44 is Ej in units of 10^^ ergs, and 10^^ci22 cm is the distance to the source. 
The duration of the neutrino burst should be of the order of the blob's ejection 
time. It should precede the associated radio outburst, or the emergence of a new 
superluminal component, that originate from larger scales 10^^ cm), by several 
hours. For sources directed along our sight line, F^^iJ'^ ~ SF^ may exceed 100. 
Thus, if the fraction exceeds a few percent, several neutrinos can be detected 
during a typical outburst from a source at a distance of a few kpc. The typical 
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angular resolution of the planned neutrino telescopes at TeV energies (e.g. Ref. I124p 
should he 6 ^ \ degree. The atmospheric neutrino background flux is '^u.bkg ~ 
10~^e~^'^j,y/cm^ssr, implying a number of detected background events A^fcfcg '-^ 
3 X lO^^(0/deg)^tdaykm~^ per angular resolution element over a burst duration 
liday day. The neutrino signals above 10 TeV from a typical micro-quasar outburst 
should therefore be easily detected above the background. 

As a particular example consider the 1994 March 19 outburst observed in GRS 
1915+105 1221, Conservative estimates of the total energy released yield -E44 = 2 — 
^44, 3^^ From the proper motions measured with the VLA, a speed of 0.92c and 
angle to the line of sight of ~ 70° arc inferred for the ejecta, corresponding to 
r ~ 2.5 and a Doppler factor 5 ~ 0.58. For these values Eq. pT|) yields a dissipation 
radius of fdiss — 10. From Fig Owe deduce that at fdiss the photopion depth t^^" 
exceeds unity at energies above 10^^ eV, where ^sy„/A = 0.1 has been adopted. For a 
maximum proton energy of 10^^ eV we then obtain ~ 0.2, and ~ 0.01(^^/0.1) 
for a 1 km^ detector. At >10 TeV this is comparable to the background, but would 
require an average of ~ 100 outbursts for detection. However, if a similar event were 
to occur in a source that is located closer to Earth, at a distance of say 3 kpc as 
in the case of the superluminal source GRO J1655-40 (which has a similar Doppler 
factor), or from a jet oriented at a smaller angle, then a few or even a single outburst 
may produce several niuon events in such a detector. Estimates of VHE neutrino 
fluxes for identified microquasars are listed in Ref. 11181 

6.4. GRBs 

As already mentioned in § 12. 3i the standard view (see Ref. [51] and references 
therein), at least until recently, was that the prompt GRB emission is produced 
behind internal shocks that form in the coasting region of a baryon loaded fire- 
ball. This scenario was originally motivated by the detection, in many sources, 
of a power law component of the prompt emission spectrum extending up to 
energies well above the MeV peak, in conflict with the ther mal spectrum ex- 
pected from an adiabatically expanding, baryon free fireball^. The radius above 
which internal shocks can form is roughly fdiss ^ 10^(roo/300)^, where Too is the 
Lorentz factor at the coasting region (see §3.1|) , and the Thomson depth is roughly 
TT ^ n'f^(TTr/Too 10^-^Cji2T^^30Q^Zlf-~^, where Eq. ^ has been used. Thus, 

'--^ 1 at r ^ a few times rdiss- Consequently, the prompt emission in this model 
originates from radii r ^ 10^^ cm. 

The difficulty with purely leptonic fireball models mentioned above may be alle- 
viated if dissipation occurs during the outflow acceleration. The claim that internal 
shocks must form in the coasting regime, since different shells cannot catch up at 
radii where the fireball is accelerating, applies only to conical geometries. Oblique- 
ness effects, inhomogeneous injection of pairs, and collision of the pair blobs with 
the surrounding matter may all lead to formation of shocks at smaller radii, where 
the optical depth exceeds unity. Bulk Comptonization may in this case lead to a 
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nonthermal extension of the spectrum weU ab ove the thermal peak. It is quite likely 
that these shocks pass through a photospherJ^HHSl^ jj^ which case the nonthermal 
photons may escape the system before being thermalized. 

The post-SWIFT discoveries of a shallow afterglow phase at early times, and 
fastly rising, large amplitude x-ray flares in the early afterglow phase " introduces 
new puzzles. According to some interpretations (e.g., Ref. I127P these observations 
imply prolonged activity of the central engine. If true, then this may indicate that 
7-rays are emitted during the prompt phase with very high efficiency, in the sense 
that the remaining bulk energy is a small fraction of the total blast wave energy 
inferred from the afterglow emission. Such episodes can be most naturally explained 
as resulting from dissipation in a pure electron-positron plasma. However, the pos- 
sibility that the shallow phase results fro m va riable micro-physical parameters and 
does not require extremely high efBciencjff^HI cannot be ruled out. 

If indeed a considerable fraction of the bulk energy is carried by baryons in the 
dissipation region, as envisioned in the "standard" model, then UHECRs and VHE 
neutrinos may be produced with a fluence comparable to the 7-ray fluence. For a 
single GRB explosion Eq. (|67p yields an integrated number of 

iV.-2W.(4)'|^WW,, (72) 

muon events, where Ej — lO^^Ej^i ergs is the total energy release in the explosion. 
The neutrinos should arrive within a time window of tens of seconds from the 
beginning of the explosion. Typically / ~ 0.1 for long GRBs that exhibit extension of 
the prompt spectrum well above the MeV peak, and m ay b e larger for bursts which 
are optically thick to pair production at these energies^. With (,p ^ /tt ^ 0.1 and 
r ~ 300, Eq. (|72p gives about 0.1 events for a typical long GRB at a redshift z = 1. 
Consequently, only nearby GRBs can be individually detected by the upcoming 
experiments. The cumulative flux produced by the entire GRB population has been 
estimated in Ref.[63]under the assumption that cosmological GEBs are the sources of 
the observed UHECRs. This requires that the rate of energy production of UHECRs 
in the energy interval 10^^ — 10^^ eV is comparable to the rate of 7-ray production by 
GRBs in the BATSE band, E - lO^''-^ erg Mpc'^ yr'^ ^2 rp^^ neutrino flux thereby 
computed is comparable to the background. However, owing to the correlation in 
time and direction of the GRB neutrinos with the associated 7-ray sources, the 
signal should be easily detected above the background, at a rate of about 10 — 100 
correlated events per year. 
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